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Abstract 

We derive formulae for some ratios of the Macdonald functions, which 
are simpler and easier to treat than known formulae. The result gives 
two applications in probability theory. One is the formula for the Levy 
measure of the distribution of the first hitting time of a Bessel process 
and the other is an explicit form for the expected volume of the Wiener 
sausage for an even dimensional Brownian motion. Moreover, the result 
enables us to write down the algebraic equations whose roots are the zeros 
of Macdonald functions. 

1. Introduction 

The (modified) Bessel functions appear in various kinds of situations. In proba- 
bility theory, for example, the modified Bessel functions of the first kind, denoed 
by Ip, appear in the explicit form for the transition probability densities of the 
Bessel processes. In this article we are concerned with the ratio of the modified 
Bessel functions. It is known that such functions represent the Laplace trans- 
forms of the first hitting times of the Bessel processes (cf. [2], [T2]) and of the 
expectations of the Wiener sausage (cf. [5]). From an analytical point of view 
Ismail et al [9l [10] have studied when such functions are completely monotone. 

We are mainly concerned with the ratios of the modified Bessel function of 
the second kind K^, so-called the Macdonald functions. Ismail [H] has shown that 
{Ky^i/ Ky){y/z) is completely monotone by expressing it as a Stieltjes transform 
of some function. From his expression we can invert the Laplace transform, but 
the resulting formula seems to be complicated. The purpose of this article is to 
rewrite the resulting ratio in a simpler form by means of the zeros of K^, and to 
invert the Laplace transform, which completes a partial result in [22]. The result 
is applied to two questions in probability theory and a study on the zeros of K^. 

Recently, in connection with the first hitting times of Bessel processes, the 
authors [7] have studied another type of the ratios of the Macdonald functions 
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and decomposed it into a sum of several functions which are easy to treat. A 
similar method via some contour integrals is effective in this article. 

The purpose of [7] is to show an explicit form of the distribution function 
for the first hitting time rjf^ to b of the Bessel process with index i/ starting 

at a. The density for r^*^^ and its asymptotics have been discussed in [8]. The 
infinite divisibility of the distribution was first investigated by Kent [12]. To be 
accurate, the conditional distribution of rj^jj under the condition that it is finite 
is infinitely divisible (cf. [ID])- General theory in the infinite divisibility of the 
distributions of the first hitting times of one- dimensional diffusion processes is 
given by Yamazato [21] . 

As referred in [7], the function K^j^i/Ky appear when we give an expression 
for the Levy measure. We may apply our result on the ratio of the Macdonald 
functions to obtain an explicit expression for the Levy measure of the distribution 

Moreover, the function {Kdi2/ Kd/2-i){fV^) in A > represents the Laplace 
transform of the expectation of the Wiener sausage for the dimensional Brown- 
ian motion associated with a close ball with radius r (cf. [5]). In the case when d 
is odd, Hamana [6] divided the function into the sum of several functions of which 
the inverse Laplace transforms can be obtained easily and deduced an exact form 
of the mean volume of the Wiener sausage by means of zeros of K(i/2-i- By using 
our result we can show that, also in the even dimensional case, the expectation is 
represented in a similar form. We should remark that the Wiener sausage for a 
Brownian motion associated with a general compact set is investigated in [T5| [T9] 
and so on, and that the same problem for a stable sausage is discussed in [Tj [T7]. 

In the results mentioned so far we have express several quantities by using the 
zeros of K^. If we consider the asymptotic behavior of the ratio {Ky^i/ Ku){w) 
as w — )■ oo or w — !■ 0, we may conversely obtain information on the zeros of K^. 
In particular, we see that the zeros are the roots of some algebraic equations with 
real coefficients and show the way to obtain the coefficients. The result is the 
improvement of that on Kn for an integer n given in [22] . 

This article is organized as follows. Section 2 is devoted to a decomposition 
of the functions K^j^i/Ky. Section 3 is devoted to a representation of the Levy 
measure of the first hitting time of the Bessel process. We calculate the expected 
volume of the Wiener sausage for the even dimensional Brownian motion in Sec- 
tion 4 and discuss its large time asymptotics in Section 5. In the final Section 6 
we study the complex zeros of the Macdonald functions. 



2. Ratios of Macdonald functions 



For each complex number v the modified Bessel function of order p is the funda- 
mental solutions of the modified Bessel differential equation 

The standard notation I^, and K^, are used to denote the functions, which are 
called the first kind and the second kind, respectively. See [21 [IH [20] et al. 
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Especially, Ky is also called the Macdonald function of order u. In this article we 
treat only the case when the all orders of modified Bessel functions are real. 

Before giving the result, we recall several facts on the zeros of the Macdonald 
function. For z/ G M let N{i') be the number of zeros of K^. It is known that 
N{v) is equal to \v\ — 1/2 if — 1/2 is an integer and that N{v) is the even 
number closest to \v\ — 1/2 otherwise. Niy) = if \v\ < 3/2. If \v\ = 2n + 3/2 
for some integer n, z^'^^e^Ky{z) is a polynomial of degree 2?7, + 1 and has a real 
(negative) zero. Otherwise, does not have real zeros. Each zero, if exists, lies 
in the half plain {z e C ; Re(2;) < 0}, denoted by C~. When N{v) ^ 1, we write 
Zy,i,Zu^2^ ■ ■ ■ ,z^,N{u) for the zeros. Since K^, is one of the fundamental solutions 
of the second order equation (12. ip . all zeros of K^, are of multiplicity one by the 
uniqueness of the solution of ordinary differential equations. This means that all 
zeros of are distinct. For details, see [201 pp-511-513]. 

Let D = {ze C\{0}; \ &rgz\ < vr} and = {z e D ; Ky{z) ^ 0} for G M. 
The purpose of this section is to show the following theorem. 

Theorem 2.1. Let w G and = max{z/, 0}. In addition, we put, for fi ^ 0, 

G^{x) = K^{xf + 7r2j^(x)2 + 27rsin(7r/i)J^^(x)J^(x), x > 0. (2.2) 

(1) When u — 1/2 is an integer, we have that, if \u\ = 1/2, 

K,+^{w) . 2u+ 



1 + (2.3) 



and that, if |z/| ^ 3/2, 



1 



KJw) w ^ Zy — w 



N(y) 



(2) When v — 1/2 is not an integer, we have that, if \v\ < 3/2, 

1 + + cos(7rz/) / — (2.5) 



Ky{w) w Jq x{x + w)G\y\{x) 

and that, if {i^l > 3/2, 

Ky+i{w) 2iy+ ^ 1 , , p dx 



, . ^ , + V + cos(7rz/) / — ,^ , (2.6) 

Ky[w) w ~^Zu,j-w Jq x[x + w)G\^\[x) 



We note that this theorem has been established in [22] when v is an integer. 
For a proof of Theorem 12. 11 it is sufficient to consider the case of ^ because 



of the formula i^T^ = and the recurrence relation 



K,^,{z)-K,_,{z) = '^-^K,{z) (2.7) 
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(cf. [201 P-79]). In fact, we see for i/ < 

K^{z) K\y\{z) z 
Formula (I2.3p is easily obtained from 

vr _,z + l 




2[z] = \l—e . 

2z z 

Moreover (12.41) has been already established in [6]. Therefore we concentrate on 
the case when z/ — 1/2 is not an integer. 

To prove fl2.5p and fl2.5p . we need three lemmas. One is an uniform estimate 
of the Macdonald function, which has been proved in [7]. 

Lemma 2.2. Let 6 G (0,37r/2) be given. For fi ^ we have 

K,{z) = ^^e-^{l + E,{z)} (2.9) 

if I arg2;| ^ 37r/2 — 5. Here \E^{z)\ ^ ^^/l-^l for a constant which is indepen- 
dent of z. 

The other two give asymptotic behavior on the real line of the functions 
involving the modified Bessel functions. Both of them are easily shown by the 
formula 

for /i ^ (cf. [201 P-203]) and we omit the detailed proofs. 
Lemma 2.3. Let C^V^^ = 0- follows that, as x oo, 

G,(x) = ^{1 + 0(1)}, ^^^^^ = e--{l + o(l)}. 



^{1 + 0(1)}, ^##^ = ^{1 + 0(1)}. 



Lemma 2.4. Let fi ^ 0. It follows that, as x oo, 

/„W-WW = -i=e-{?^iii + o(i)^ (2^11) 



We are now ready to show Theorem 12.11 We only consider the case when 
z/ — 1/2 is not an integer. Let a G (0, 1) and w G D^. For z G with z ^ w, we 
set 

1 K+,{z) 



z°'{z-w) Ky{z) 
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Letting e and R be positive numbers with e <1 and e < R and setting 



7r 



On^e = Arcsin — G (^0, - 

we consider the same piecewise C^-curve 7 as in [7] defined by 

70 : ^ = Re'', -71 + eR,e ^e^ix- Or,,, 

71 : z = X + ie, —RcosOr^ ^ x ^ 

72 : 2; = X — ze, —RcosOr st ^ x ^ 

73 : 2 = ee^^ -7r/2 ^ ^ ^ 7r/2, 

7 = 7o + 7i - 73 - 72- 

We take R so large and e so small that w and all zeros of K,^ are inside 7. Then, 
setting 

for = 0, 1, 2 and 



27rz 



73 



we have 

^{R, a, e) = %{R, a, e) + !P'i(-R, a, e) - !Z^2(^, a, - 5^3(a, 



(2.12) 



The singular points of /^^^ inside 7 are i/; and the zeros of K^, which are all 
poles of order one. Hence the residue theorem yields 



Res(^/;;/- ) 



iiN{u) = 0, 



W{R,a,e) 



N{u) 



Resiw; /-J + Res(z.,,; /-J if Niu) ^ 1. 



(2.13) 



Here Res(f ; /) is the residue of a function / at a pole v. It is obvious that 

1 K,+i{w) 



Res(^; /- ) 



w" KJw) 



When N{y) ^ 1, by the formula zKl{z) — vKy{z) 
we have 



Res(2jyj', fyn) 



1 



Ky+l{Zi,^j] 



-zK,+i{z) (cf. [201 P-29]), 
1 



Hence we obtain from (12.131) 



lim lim lim \I/{R,a,e) = < 

e^O a— >0 _R— s>oo 



KJw) 



N{u) 



if iV(z/) = 0, 
if N{v) ^ 1 



(2.14) 



if we show that the hmit on the left hand side exists. 

We fix e > 0, a > and consider the right hand side of (12.121) . By (12. 9p we 
have for sufficiently large R 



R'^R-lwl I-AJR 
which immediately yields that 

lim %{R,a,e) = 0. 
For the integral \Pi{R,a,e), we have 

^i(R,a,e) = / 7 r ^7/ r^dx 

1 f^'°'^R>^ 1 K,+^{-x + ie) 



2'Ki Jq {—x + ie)°'{—x + ie — w) Ky[—x + ie) 



dx. 



Note that 

for z G D and /i ^ (cf. [201 p.80]). Then, setting 

we have 
and 

?Z^i(it, a, e) = / — — 7 rcta; 

^ 27r« 7o e*™(x - - «5 + w) 

' ^ '-'^'^ -dz. 



(2.16) 



27rz J^o e*™z"(z + u;) 



where 7° is the line in D defined by 7° : z = x — ie, ^ x ^ RcosOn^e- 
We now define three paths as follows: 

'yl : z = ee'\ -tx/2 ^0^0, 

7^ : z = X, e ^ X ^ R, 

ll: z = Re'', -9n^, ^9^0. 

Since w is inside 7, we have that | lm{w) | > e if Re(ty) < 0. There is no zero of K^, 
on the real axis and the integrand of the right hand side of (I2.16P is holomorphic 
inside and on the contour consisting of 75", 7J, 7^ and ^f. Then it follows from 
the Cauchy integral theorem 

^i{R, a, e) = ^l{a, e) + ^^{R, a, e) - ^^{R, a, e), 
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where 

fO 



Wt{a,e) = — / , , — -dd, 



' ' ' 2711 e*™x°(x + w) 



If -7r/2 ^ ^ ^ 0, we deduce from ([23]) and fl27[5|l 



for sufficiently large R and 



K^{Re(<>+^)) 



< 



1 + 
1 - 



which immediately implies 

lim ^^{R,a,e) = 0. 

R^oo 

Moreover, from fl2.17p for R = x and ^ = 0, it follows that 

lim ^^{R,a,e) = — H —^^4^^L_dx, 

for which we write 

In order to consider \p2{R, we recall 

K^{e-'^z) = e'-''K^{z)+t'Kl^{z) 

for z ^ D and /i ^ (cf. [201 P-80]). Then we have 



where 



-^i.,2(2), 



In the same way as ^i{R, a, e), we can show that 



where 



1 /"^/^ £e^^(7,,2(£e^'l 



^l{R,a,e) = — r ^3u^ , 



^^(R,a,e) = — , ^'^'.^ ^ — -de. 
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(2.17) 



Similarly to (-R, «, e) and \l/f{R., a, e), it is easy to see that 



lim iP'a (-R, a, e) = IP's (a, e), lim (Z^s (-R, a, e) = 0, 



where 



1 r ^?-,2(x) 



■2 ' 



27ri e-*™x"(x + w) 



Therefore we conclude that 

lim ^{R,a,e) = ^l{a,e) + ^i{a,e) + (a, e) - - ^3(0, e). (2.18) 



We next consider the limiting behavior of the integral on (12.181) as — ?■ 0. 
We first let a — )■ 0. Then it is easy to see 

Iim5^i(a,e) = — / jTT, de, 



a^O 27T Jq 



ee''^ + w 



\imW2,{a,e) = — / — — - — ^— db'. 



It is known that 



log-Vl + o(l)} if/i = 0, 

iA/^/- {1 + 0(1)} if/x>0 



2 V^. 

as \z\ — in D (cf. [201 P-512]). Moreover, it is easy to see 

1 + 0(1) if/i = 0, 



r(/i + i) V2 

as \z\ — )■ in D by the series expression. With the help of (I2.19P and (12. 20 p . we 
obtain that zgy^i{z)/{z + w) and zgy^2{,z) /{z + w) tend to if z/ = and to 2u/w 
if z/ > as |2;| — >■ in D. Hence we have 

lim lim ^}{a, e) = lim lim \PUa, e) = — . (2.21) 

Moreover, since zKyj^i{z) / {z — w)K^{z) converges to —2v/w as \z\ — > in D, we 
obtain ^ 

lim lim (^'^(a, e) = . 

and 

lim \imW}ia, e) + !P'?(a, e) - ^^{a, e)} = —. (2.22) 

e-5-0a-s>0 W 
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Net we set 



Then, putting 



and 



+ t7Te-'^''{e''"'h^i{x)K,{x) + e-'""K,+i(x)/,(x)} 
+ t7re'^''{e-'^''h+,{x)K,{x) + e'^''K,^i{x)h{x)}. 



F^{x) gu,i{^) 9v,2{x) 

K[x) 



Gy{x) e 
where Gy{x) is given by (12.21) . we have 



Recall the formula ^ 

Ki,+ l{x)ly{x) + IyJ^l{x)Ki,{x) = - 

X 



(cf. [201 p.80]). Then we get 



e-^™{e*'^^J,+i(x)JC(x) + e"'^''K+,ix)Ux)} 

= e-^™<^ + 2ism(7iiy)I^+i(x)KJx) 

I X 

+ e^™<^ 2ism(7iu)I^+i(x)KJx) 

{ X 

2cos7r(a + z/) ^ . , , . , ,^ / n / n 
= |-4sin(7ri^) sm[7ra)I,y+i[X)Kiy[X). 

X 

Hence, letting 

= 2isin{7ia)K^+i{x)K^{x), F^.2(x) = 2z7r^ sin(7ra)/j,+i(x)/^(a;), 
2z7r cos Trfo^ ~h ^) 

F'u,3{x) = , Fi,^4{x) = 4i7rsin(7rz/) sm{7ia)Iu+i{x)Ki,{x), 



we have 



We set 



1 r 1 



e) = — ^ . ;^'7 / cia:, 1 ^ A; ^ 4. 



By virtue of Lemma 12. 3[ the functions 

1 K,+^{x)K,{x) 1 



x + w Gy{x) ' X + If G',/(a;) 
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are integrable on (e, oo). Hence we get 



lim^i(a,e) = lim^4(a,£) = 0. (2.23) 

a— 5>0 o— !>0 



The integral ^3(0, e) is written as 



+ w)G„{x) 



and, with the help of Lemma 12.31 we can derive 

x{x + w)G^{x) 



It follows from f l219|) and fl2:20|l that 



dx 

lim^3(a,e) = cos(7rz/) / — (2.24) 



log- {l + o(l)} ifz/ = 0, 
G'.(x) = <^ V (2.25) 

'^{1 + 0(1)} if.>0 

as X — 7- 0, where = l/4'^~^{r'(z/)}^. This implies the convergence of the right 
hand side of (12.241) as e — j- and 

f°° dx 

hm lim <p3(a, e) = cos(7ri^) / — , ^ , , . (2.26) 

e^Oa^O ' ^ Jo x{x + w)G^{x) ^ ' 

It remains to consider ^2 («,£), 

<f.( \ ■( N r 1 /.+i(3:)/.(x) , 

<P2(a, = TT sm(7raj / — ^ , , dx. 

We should remark that the function 

1 /^+i(x)/^(x) 
X + U7 Gy{x) 

is not integrable on (e, 00). We write 

^2(a, e) = -^2(«. ^) - ^2(". ^) + ^2(«. ^) + ^2(«. 

where 



, sin(7ra) p 1 i^.(x)2 

^^(a,e) = — ^ — -\ — r ^ , , dx, 

^ x°(x + u') G^(x) 

^2/ N / N / N P 1 /^+l(x)i^Jx) , 

(Pnla, e) = 2 sm(7ri^) sm(7ra) / — ^ , , ctx, 

^ ^ ' Js x'^ix + w) Gy{x) 

^3/ ^ . . ^ r 1 /.(x){4(x)-4+i(x)} ^ 

<Pn a, £ = TT sm vra / — „ , , dx. 



2V"^' ' 



x"(x + w) G,y{x) 



^A, , sin(7ra) /"°° dx 
^^(a,e) = 



TT 



x"(x + w) 
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We can easily deduce from Lemma [2.31 

lim $l{a, e) = lim ^l{a, e) = 

in the same way as fl2.23p . 

To compute ^l{a,e), we note the following: 

/.(x){J.M - J„..(x)} ^ 2^{i ^ ^ ^ ^_ 

which is obtained from (12.101) . (12. lip and Lemma [2. 3 [ Since 

1 h{x){h{x) - h+i{x)} 
x + w G^{x) 

is integrable on (e, oo), we have that 'P\{a,e) tends to as a — )■ 0. 
The calculation of ^2(0^)^) is easy. In fact we have 

^4, , sin(7ra) /"°° dx wsmina) f°° dx 



sin(7ra) wsin(7ra) dx 



1, a 



for any £ > 0. Hence we conclude from (I2.23P and (I2.26P 

limlim(p(a,£:) = 1 + cos(7rz/) / — , (2.27) 

e^oa^o ^ ' ^ io x{x + w)G^{x) ^ ' 

Combining (EUD, ([218]), ([221]), fl2:22|l and ([227]), we complete our proof of 
fl^ and ([22]) in the case of ^ 0. 



It may be worthwhile to mention that the method used to prove Theorem 12.11 
can be applied to the decomposition of Ky^pjKy, and we can show the following. 

Theorem 2.5. Let v and p be real numbers with z/ ^ 0, —u ^ p < 1, p ^ 0. 

In the case when there is no integer n such that z/ = 2n + 3/2, we have that, if 
V < 3/2, 



Ku{w) Jq x + w Gy{x) 

and that, if v > 3/2, 



dx 



^ z^j - w K^+i{z^j) Jq x + w G^{x) 



where 



Hu,p{.x) = - cos7r(z/ + p)Ki,+p{x)I^{x) 

+ cos{7Tu)I„+p{x)Ki,{x) + ^^^^^P^ K^+p{x)K„{x). 
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When z/ = 2ra + 3/2 for some integer n, we have that 



N{u) 

1 + A^,p + ^ 




where 



A 



lim 

5-5.0 



sin(7rp) f 1 Ky+p{x) 



dx, 




^ J\x-xo\>5X + wK^{x) -7ll^{x) 



ifzeD, 



if z E {—oo, 0) 



K^+i{-z) + Trlu+i{-z) 



and Xq is the unique real zero of Ky. 

With the help of Kp = K_p and fl2.8|] , we furthermore see that Theorem 12.51 
gives a representation for Kp/K^ for every /i, z/ G M. 



For z/ G M the one- dimensional diffusion process with infinitesimal generator 



is called the Bessel process with index u. The classification of boundary points 
gives the following information. The endpoint oo is a natural boundary for any 
1/ G M. For z/ ^ 0, is an entrance and not exit boundary. For — 1 < < 0, 
is a regular boundary, which is instantly reflecting. For z/ ^ —1, is an exit but 
not entrance boundary. For more details, see [TT] and [IB] for example. 

For a, 6 G M let rj^jj be the first hitting time to b of the Bessel process with 

index z/ starting at a. The conditional distribution of rf*^^ under t^'^^ < oo is 
infinitely divisible. The purpose of this section is to give the exact form of the 
Levy measure m^al when ^ 6 < a by applying Theorem 12.11 

It is known that, when ^ a < 6, the distribution of r^^'J is a mixture of 
exponential distributions. Let us recall the results in [2]. See also [H]. In this 
case, the Laplace transforms of the conditional distributions are given by the 
following. For A > 0, if 6 > and z/ > —1, 



3. The first hitting time of the Bessel process 




E[e 



-At 



0,b 




if < a ^ 6 and z/ > — 1, 
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if < a < 6 and z/ < -1 



hj /_,(6V2A) 



Combining these results with the formula 



I, Ax) 



X \ 1 -i-r / X 



for ii> —1 and x > 0, where {j^,n}^i is an increasing sequence of positive zeros 
of the Bessel function of the first kind of order /i, we obtain the following 
expressions for the Levy measures: if 6 > and v > —1, 



' n=l ^ 



dx X ^-^ 

n=l 



— ^ ' ^ e 



if < a ^ 6 and z/ > -1 

^"^SCfel^) _ 1(0,00) (a;) 



X 

n=l 



if < a < 6 and u ^ -1 

^a,bi^> i(0,oo)W ^/ _i 



^^^il(^) 1 (0,00) (X) ^ 
.7T 

n=l 



X 



I e 262 _ g 2a 



respectively, where 1a is the indicator function of a set A. 
The following is the main result in this section. 

Theorem 3.1. For ^ b < a the support of the Levy measure m!"^\ is [0, cxo) 
and it is absolutely continuous with respect to the Lebesgue measure. We have the 
following expressions for the density p^l, x > ofm^^l. 
(l)//a>0, 



(2) If a > 0, u - 1/2 and u ^ -3/2, 



Mr \ 



N{u) 



V^irx^ 2y/7rx^ ~[ Jo 
(3) lfa>0, -3/2 <u <0 andu^ -1/2, 



cosinu) r r 1 



V2nx^ 2Vnx^ Jo Jo vG\u\{v) 
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(4) If a >0, u- 1/2^ Z and u < -3/2, 



V^TTx'^ 2^A^a? ~^ Jo 



COS TTZ/ 




2\/ TTX^ Jo Jo 



(5) IfO<b<a andu = ±1/2, 



a — b 



V2' 



(Q) If < b < a, u - 1/2 e Z and \u\ ^ 3/2 

Niu) 



\/2nx^ 2V nx^ ~^ Jo 

(7) //O < 6 < a, ^ \u\ < 3/2 an(i ±1/2, 
a-b cosfTTz/) /•°°/'°° 1 



/ e 4^ ( e — e v^*" j (i^. 



If Q <b < a, V -1/2 ^Z and \v\ > 3/2, 



Pali^) 



a — b 



N{v) 



y/27rx^ 2\/ irx^ ~[ Jo 



± 



cos TTZ/ 



1 _ti ( „_S2L __S2L \ 

, , — — -— e 4x e _ e id^drj. 

2V7rx3 Jo Jo ^<^Ih(^) V / 




The rest of this section is devoted to the proof of Theorem l3.1[ By the formulae 
for Laplace transforms for rj^), we have, for A > 0, if a > and u < 0, 



kio < ool 



r{\u\){aV2Xr 



K^{aV2X) 



if < 6 < a and u G 



E[e -l4Uoo]=^-J — 



:|.|(&V2A) 

We represent K^{x) by means of and zeros of Ky. 

Proposition 3.2. For x > we have the following formulae. 
(1) 7/^ = 1/2, 

1 TT 

log{x^/^i^'l/2(x)} = - log - - X. 



(3.1) 



(3.2) 



(3.3) 
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(2) ///i - 1/2 e Z and fi ^ 3/2, 



log{x'^J^^(x)} = log{2'^-ir(/i)} - X - ^ log (3.4) 

j=l ^A'-J' ^ 

(3) //O < /i < 3/2 and/i 7^ 1/2, 

logKi^^(x)} = log{2^-ir(/i)}-x-cos(7r/i) / — — -log^ dy. (3.5) 

(4) If fi - 1/2^ Z and > 3/2, 



log{x'^i^^(x)} =log{2'^-ir(/i)} - X - ^ lo, 



3=1 ^ (3.6) 

-cos(7r/i) / ——log rfy. 

Jo yGf^iy) y 

Proof. Formula (13. 3p is obtained from the explicit expression for Ki/2{x)- 
In order to show the others, we note the following formula: 

— log{x^i^,(x)} = - + — , (3.7 

dx J^u\x) X 

which can be derived by ^1} and [x^'K^^ix))' = -x^K^„i(x) (cf. [201 P-79]). 
If z/ - 1/2 G Z and z/ ^ 3/2, it follows from (EID and ([32D that 

d ^^""^ 1 

-log{a:^K.(x)} = -l-5:-— 



■ 1 "j^i? 



Then we obtain that, for any £ > 0, 



log{x'^if,(x)} - log{£^J^,(£)} = -[x-s) + Y. f TZ: 

7 = 1 ^ 



Hence, letting e — 0, we get (13. 4p with the help of 

lim x^K^i.x) = 2''-^r{v). (3.8) 

If < < 3/2 and z/ 7^ 1/2, it follows from and (jSTD that 
log{x'^-ft'j,(x)} = — 1 — cos(7rz/) ^ '^^ 



dx " Jo y{y + x)G^{y)' 

Hence we have 

\og{x''K^{x)} - \og{e''K^{e)} = ~{x - e) - cos(7rz/) ^ ^ 



dy 



y{y + OGu{yy 
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Note that Gi,{x) is positive for x > unless v is 2n + 3/2 for any integer n. Thus 
it follows from fl3.8D that 



logix"" K^{x)} - log{2^-^r(i^)} = -x- cos(7rz/) / 



X POO 



dy 



'0 ^0 

By the Fubini theorem, the right hand side is equal to 

r dy r di 

-X — cosyKu) 



yGu{.y) Jo y + C 

which yields (13. 5p . 

We can show fl3.6p in the same way. □ 



In order to see Theorem 13.11 we only have to check 

02(A) := \ogE[e-'<^ \ rij < oo] = / (e"^^ - l)pi'l{x)dx 

Jo 

for each case. It follows from (13. ip and (13. 2p that, if a > and u < 0, 

0l:J(A) = log{(av^)l'^lir|.|(av^)} - log{2l^l-ir(|H)}, 
and that, if < 6 < a and G R, 

0j;)(A) = log{(av^)l'^lK|.|(ay2A)} - \og{{bV2XpK^,^ibV2X)}. (3.9) 
The following lemma gives rise to Theorem 13.11 for u ^ hj Proposition 13.21 
Lemma 3.3. Let c > 0, u > and z G C~. For X > it follows that 

2A = - / dx, (3.10) 



3 



and 



v2vrr 

log j== , r^ dx I e~^^vr.d^ (3.11) 

z — cv 2A Jo 2v TTX'^ 



7] + cV2\ 1 
V vGuiv) 

oc^-\x_i roc POO ^ _e_AlL 

d^ / „ , e *^ -y^cd^di]. 



(3.12) 



2 a/ ttx^ Jo Jo 

Proof. We recall the following formula (cf. [3i p. 361]): for — 1 < p < 0, 

POO 

/ y'-\e-y-l)dy = r{p). 
Jo 

Setting p = -1/2 and noting r(-l/2) = -20^, we obtain fIXTU]) . 
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Recall the formulae 



roo 

e~—^x-'/'dx = ^e-^^\ (3.13) 



'0 



e--^x-'"^dx = 2 / e^^di/ = A, ^ > 0. (3.14) 

Jo 



and, for G C with Re(a) > and Re(/3) > (cf. [31 p.361]) 



Then we obtain 




JQ 



oo -ax _ p-l3x o 

dx = log ^ (3.15) 

X a 



——e-t^-'^^dxd^ = 2v^ / d^ 

^ Jo ^ 

VX + a 



2v^log. 



a 



Hence, applying the Fubini theorem to the right hand side of (13. lip , we obtain 

1 ^ .. e"^^ - 1 «2 . ^ e~^^ - 1 



e^^d^ I — - — e i^dx= / e^<= d^ 



log 



The calculation of the right hand side of (I3.12p is similar. Since Gy[x) is 
positive for x > 0, we may apply the Fubini theorem and, applying (I3.13p . (I3.14p 
and (I3.15P again, we get 



/ -^^dx / / ^ , , e ^'^d^df] 

Jo 2Vnx^ Jo Jo vGuiv) 

dr] _jiL. ,^ /""^ e-^^ - 1 _£i , 

e ^'^dt, — e 



' ' e dC, 



vGuiv) Jo ^ 

r 1 , V/V2c 

log —= -j^dr]. 



Jo r]G^{7]) y/X - ri/^/2c 
This immediately implies (I3.12p . □ 

Remark 3.4. If = 0, it follows from flT^ that the left hand side of flXT^ 
diverges. 

We finally consider the case of = 0. Since Kq{x) = —Ki{x), we have 

d K,{x) r 

\ogKo{x) = --—- = -1 



dx Kq{x) Jq ri{r] + x)Go{r]) 
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by Theorem 12 ■![ Hence, by (13. 9p . we get 

i>a\/2X POO J 

<P^^l{X) = -{a - b)V2X - ' 



bV2\ Jo 

r]{r] + x)Go{r]) 

-{a - b)V2\ - / ——-log -=dr]. 

Jo vGoiv) r] + bV2\ 



(3.16) 



Note that 

1 T] + ay/2X 

vGo{r]) 7/ + bV2\ 
is integrable on (0, oo) by (12.251) . For c > and e > 0, we write 



' 7] + cV2X 1 , /""^e-^^-l dr] 

^ vGoiV) Jo 2Vnx^ Je vGoiv) Jo 

This formula immediately implies 



r 1 , r/ + aV2A^ 
/ ^ / N log i=dr] 

Je vGoiv) V + bV2X .3 



Jo 2v^ Je vGo{v)Jo V / 

If we show that the integrand is integrable on (0, 00) x (0, 00) x (0, 00), we can 
conclude 



£^0 



V + cV2\ 1 ^ 
hni / log TTT^dr] 



= / ; dx I ^ , , / e 4^ e ^6 _ e ^« 

Jo 2v^ Jo vGoiv)Jo V y 

and obtain Theorem 13.11 in the case of z/ = 0. 

Since A > and < 6 < a, the integrand on (I3.17p is non-negative. We have 



that 



^ p e V2b — e 

2Vvrx3 r/G'o(r/) V 

1 — e~^^ 1 _ii __S22_ / __SlL + _S22_ 

= g 4a: g j X — 6 

2v/^ ^^0(7/) V 
1 _g-A^ ^ 1 

\ g 4x g v^a I — — — 

~ 2V27rx3 Gol^) 



From (I3.12P and (I3.13P we obtain 

fOO pCO n c„ 1 ^— Ax 




Jo 



^g-ijg ^^—d^dx = 2y/^ / e ^(1 -e-«^^)(i^ 

^2v/^/ 
Jo 



SzL ,^ 2V27ra 

g V2ad^ = . 
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Hence the integral 

foo coo g— Ax 




Jo 2^fK3^ vGoiv) 



is bounded by a constant multiple of l/riGo{r]), which is integrable on (0, oo), and 
fl3.17|) converges as e — )■ 0. 

4. The expected volume of the Wiener sausage 

Let r be a given positive number. The Wiener sausage {W{t)}t^Q for the Brow- 
nian motion with radius r is defined by 

W{t) = {xeR'^; x + B{s) G U for some s G [0, t]} 

for t ^ 0, where {B(t)}t^Q is a Brownian motion on M*^ and U is the closed ball 
with center and radius r. For t > let 

L{t) = I P^[r ^ t]dx, 



where r = inf{t ^ 0; B{t) € t/} and is the probability measure of events 
related to the Brownian motion starting from x G M''. It is easy to see that the 
expectation of the volume of W{t) coincides with the sum of L{t) and the volume 
of U. 

In the case when d is odd, the explicit form of L{t) has already given. One 
and three dimensional cases are easy. Indeed, we have that 



L{t) 



2^2tJ^ ifrf=l, 
2nrt + Ar'^y/2Td if = 3. 



These formulae can be obtained directly from the well-known formula for 
Px[t ^ t]. For details, see O [11], [15]. In the higher dimensional cases, the 
authors [7j recently obtained an explicit form of Px[t ^ t]. However it is not of a 
convenient form for the integration on x. 

We here consider the Laplace transform of L given by 



e-^'L{t)dt = ^^^^^4= "^^^ ^ , A > 0, (4.1) 

i^,/2_i(rv^)' 

where S^-i is the surface area of c? — 1 dimensional unit sphere (cf. [S]). When 
d is odd, since Kd/2{x) / Kd/2-i{x) may be expressed by the ratio of polynomials 
for X > 0, and the right hand side of (14. ip may be represented by the linear 
combination of rational functions of the following four types: 



v^' A' Vas' ^ 
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Hence the Laplace transform on (14.11) can be inverted. When d is odd and more 
than or equal to five, Theorem 1.1 in [6] shows that, for t > 



L{t) = Sd-ir 



d-2 



{d - 2)t y/2r^ 
2 ^ d-4 



j=l l^j ; 



r X 

e 2t ' -j 



. (4.2) 



Here we have used Zj^^ and A^'^ instead of Zd/2-i,j and N{d/2 — 1), respectively. 

Our goal in this section is to give similar results in even dimensional cases by 
applying the results in Theorem 12.11 

Theorem 4.1. For x > we set G^^\x) = Gd/2-i{x)- 
(1) If d = 2, we have 



2t , V2r^ xy - I + e-'^y _r2£ 



Tt Jo Jo y'G(^Ky) 



L{t) = 27rr \l — + 

(2) If d = 4, we have 

L Jo Jo y^Gi^Ky) 

(3) If d^ 6 and d is even, we have 



e 2t dxdy 



V2r^ r°°r°° 1-e-^^ 



e 2t dxdy 



L{t) = Sd^ir 



d-2 



{d - 2)t V2'i 



,3 



d-4 



TTt {zf^f JO 



r x_ 

e 2' ' 



+ 



-xy 



3 poo POO 

Jo y'G^'Ky) 



e 2t dxdy 



We set T{t) = L{2rH) for t > and then deduce from that, for A > 



POO 

/ e-^'T{t)dt 
Jo 



A3 IQ/2-iiVX) 



For a proof of Theorem 14.11 we consider the function: 

/C(VA) 

Let Tjy be the inverse Laplace transform of Then we have the following, which 
immediately yields Theorem 14.11 

Theorem 4.2. Setting = max{i^, 0}, we have the following. 
(1) If \u\< 1/2, 



rr. , . + it COS TTZ/ /•°°/'°° 

^' y/lTt Jo Jo y''G\,\{y) 



^ + '''e-^cixdy. (4.3) 
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(2) //l/2<|z/|^l, 

(3) //1< li^l <3/2, 



cos TTZ/ 



e 4t dxdy. 




vrt Jo Jo y^G\v\{y) 



%{t) = 2z/+t 



1 



cos vrz/ 




-xy 



2(|z/| - 1) 7o Jo y''G\,\{y) 

(4) // |z^| > 3/2 and v — 1/2 is not an integer, 



T^{t) = 2v^t 



1 



-y- 



2{\u\ - 1) ^ JO 

cos(7rz/) /•°°/'°° e-^s' 



vrt Jo Jo y^G\u\{y) 



e 4t ^''•"'i^'dx 
dxdy. 



(4.4) 



e *t dxdy. (4.5) 



(4.6) 



From now on, we shall treat Kyj^i/K^ as the function on the half real line 
(0, oo). For our purpose we need three lemmas. 



Lemma 4.3. We have that, as x 4 0, 







r 2i^+ 

x^ 
2u+ 






Ky+l{x) 


= < 


X 






X 

2p+ 


^2(|H- 

X 


1) 






. X 


^2{\u\- 


1) 



3 

+ o{x) ifl<\v\<-, 

+ o{x') ^f\u\>^,u-^^Z. 



In order to prove Lemma 14. 3[ we need to show an asymptotic behavior of 
K^{x) as a; I 0. 

Lemma 4.4. Ifn + l/2<fi<n + 3/2 for an integer n ^1, we have 



as x J, 0. Especially, we have 



{2xf + o{x 



2n+l\ 



2 \x 



- e^^'i 1 + X + - 



12fi-3 2 12/i-5 



X + -- 



22^- 2 6 2/i-2 



-a;'' + oix" 



Proof. It is well-known that, for /i > —1/2 




2xr(/x + 1/2) 



{2xyr{^i + i/2) 7o 



e "y"^ 



2x 
2x\ 

y J 



dy 



dy 



(4.7) 
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(cf. [m p.l40], [201 p.206]). The Taylor formula yields that 
e-V^^^(^l + ^j dy 



oo 







k=l 



2n + 1 J \ y J \ y 
for some ^ G [0, 1]. Since 2fi — 2n — 2 > —1. Since fi — 2n — 3/2 < 0, we have 

\ y J 

/•oo / 9tT,\ A*-2n-3/2 



V y 

oo 



/■oo 

S / e-V^"^""^t^2/ = ^(2/i-2n-l). 
io 



Hence the dominated convergence theorem yields that 

/•oo / o \ A'-2n-3/2 

lim / e-^|/2/.-2n-2 1 + ±t rf^/ = r(2/i - 2n - 1). 

Jo \ y J 

Therefore we obtain 

r(2/i) + (2/i - l)r(2/i - l)x 



'"'^^'^^ (2a;)^ r(/i + 1/2) 



2?i ^ _ 1\ 



,2n + 1 

J, 0. With the help of the formula 



^ 2 ^ -2n- l)(2x)2"+i{l + o(l)} 



2'^^^r(z)r(^z + =v^r(2^) (4.8) 

(cf. [m p.3]), we easily obtain the assertion. □ 



Proof of Le'mma \4.S\ We only consider the case of z/ ^ 0. If we show this case, 
the result for u < follows from fl2.8p . 

If > 3/2 and — 1/2 is not an integer, it follows from fl4.7p that, as a; J, 0, 

K,{x) X ^ 2{u-l)^''^'' 
22 



When l</i<3/2, we use the formula 



e-V^"Ml + ^J dy 







oo 



\ ^ J \ y J 



in a similar way to Lemma 14.41 Then we have 



/ A f2Y 12^-3 2 , 2^ 



Combining it with (14. 7p . we deduce that, ii 1 < u < 3/2, 

+ o(a;). 



K^{x) X 2(z/-l) 

In the case of 1/2 < fi ^ 1, the calculation is simpler. Indeed, by the same 
calculation as (14. 7p . we can easily obtain 

K,{x) = ^ 0^ V^{1 + X + o{x)} (4.9) 



in virtue of the formula 



11.^1/2 



°°e-V^-^fl + -) dy 

V y 



Hence we deduce from (14. 7p and (14.91) that 

KJx) X 1 + X + oix) X 



We finish the proof of Lemma 14.31 □ 
For an integer /c ^ 1 we set 

^ dy 



Cy,k — ^ I, , Qv,k — I 

j=i Jo 



In virtue of Theorem 12.11 and Lemma 14.31 we can derive the connection between 
C,.,fc and g^^k- 
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Lemma 4.5. (1) If 1/2 < \u\ ^ 1, we have 

1 + cos(7rz/) = 0. (4.10) 

(2) //I < |z/| < 3/2, we have 

1 + g^^i cos(7rz/) = 0, 

^' ^ ' 2(|z/|-l) 

(3) // \v\ > 3/2 and u — 1/2 is not an integer, we have 



(4.11) 



1 + Ciy.l + Qiy.l COs(7rZ/) = 0, 

. N 1 (4-12) 

G,2 - Qu,2 cos(7rz/) = ^^^^^ _ 

Proof. Recall that, for z/ 7^ 

{71 
— e^''{l + o(l)} as X cx), 
^-^{1 + 0(1)} asx;0. 

If 1/2 < |z/| ^1, we obtain by fl4.13p that l/y^G\,^\{y) is asymptotically equal 
to as y I 0. This implies the convergence of g^^i. The dominated 

convergence theorem shows that, as x J, 0, 

^,,,1 + o(l). 



'0 (2/ + a;) 2/^1^(2/) 
Hence we deduce from (12. 5p that 

. s = 1 H h cos(7rz/)^^,i + o(l). 

With the help of Lemma [4. 3 j we conclude (14.101) . 

In the case of 1 < |z/| < 3/2, l/y^G\i,\{y) is integrable by (14.131) and we have 

dy [1 X ^ dy 



{y + x)yG\^\{y) Jq [y y{y + x) j yG\^\{y) 

dy (4.14) 



= Qu,i - xgu,2 + o{x) 

by the dominated convergence theorem. Then we get by (12.51) that 

K^+i{x) 2z/+ f \ , f \ 

, , = Vl^Qvx cos(7ri/) - XQy 2 cos(7rz/) + o[x) 

Ky[x) X 

and hence, combining with Lemma [4.3[ we get (14. lip . 
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If |z/| > 3/2 and z/ — 1/2 is not an integer, we deduce 

N(u) N{u) 



Niu) ^ 

= Ci +xY^ — — 

j=i ^'^jy^'^'j 

= Cu,l + xCu,2 + 0{x). 

Since l/y^G\,y\{y) is integrable on (0, oo) by fl4.13p and fl4. 141) is valid for this case. 
Therefore it follows from fl2.6l) that 

„ . . = h 1 + Cu,i + Qu,i cos(7rz/) + a;{Ci/,2 - 0u,2 cos(7rz/)| + o(x). 

We immediately obtain (14.121) by Lemma 14.31 □ 

Remark 4.6. In the case when \u\ > 3/2 and z/ — 1/2 ^ Z, we can derive 

C!^,3 + Qk3 cos(7rz/) = 0. 
It is not necessary for the proof of Theorem 14. 2[ 

In virtue of Theorem 12.11 we have that, if |z/| < 3/2 and — 1/2 ^ Z, 

= ^ + ^ + cos(7rz/) / ^ \ — (4.15) 



' 'Jo V^\y^+y)yG\,\{y) 
and that, if \u\ > 3/2 and z^ - 1/2 ^ Z, 

2z.-^ 1 1 



+ cos(7rz/) 



oo 



dy 



(4.16) 



y^{VX + y)yG\,\{y)' 



Lemma 4.7. Fort > let 



rr xy-l + e-'^y , , , , 

qu{t)= / — p{t,x)dxdy, 

Jo Jo rGlHl^) 

where 

1 

j)(i:, x) = —;=e 4t . 

Then, we have, for A > 

e-^'q,{t)dt= . (4.17) 

Jo VA^(VA + y)yG'i^i(y) 
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Proof. We first recall the elementary formula 



' 1 
e-^'p{t, x)dt = A > 0. 

Then, with the help of the Fubini theorem, we deduce 



(4.18) 



e-^^q^{t)dt 



dy 



dx. 



lo y'^G\v\{.y) 

Carrying out the elementary integral in x, we obtain fl4.17p . □ 
We now complete our proof of Theorem 14.21 If \v\ < 3/2 and — 1/2 ^ Z, 

T^{t) = 2z/+t + + cos{7ru)q^{t). (4.19) 

When |z/| < 1/2, (14.191) immediately implies (143|) . 

When \u\ > 1/2 and z/ - 1/2 ^ Z, we have by fHl3|l that l/y'^Gi^i{y) is 
integrable on (0, oo). Since 



1 - e-^^ 



y^G\u\ (y) 



pit, X 



< 



xp{t, x) 



we have 



q,{t) = 2J-g, 



y^G\^\{y) 



'''^ Jo Jo y'^G\^\{y) 




p{t, x)dxdy. 



(4.20) 



Combining this formula with (I4.10p and (I4.19p . we conclude (14. 4p . 

When > 1 and z/ — 1/2 ^ Z, we can further improve (I4.20p . Indeed, since 
l/y^G\y\{ii) is integrable on (0, oo), we get 



oo /"OO 



-xy 



io y'^G\^\{y) 



pit, x)dxdy 



converges and we get 




-xy 



Jo y^Gi^iiy) 



p{t,x)dxdy. (4.21) 



Hence, we get (143]) by (gUD, (1421]) and Lemma Ol 
To invert in the case of |z/| > 3/2, we note 



1 b 



Then it follows that 

Niu) 



a±b a a? a'^{a±b) 



N{u) 



X 



j=l - 



j=l ^i/jA^vj 
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Hence f l4.16p is equivalent to 



oo 



dy 



+ cos[txv) , 

'o VA3(VA + y)i/G|,|(y) 

With the help of (I4.18p . we easily see that, for z G C~ 

1 



/•OO /'OO 

/ e'^^dt I e'''p{t,x)dx 
Jo Jo 



Hence, we deduce (14.61) from fl4.2ip and Lemma 14.31 

5. Large time asymptotics of the Wiener sausage 

This section is devoted to show an asymptotic behavior of L{t) for large t in even 
dimensional cases. Le Gall [15] considered the Wiener sausage associated with a 
general compact set and proved 

[c[^)t + cf + cf t2-d/2 + o(ti-'^/2) if ^ 5 

and gave the explicit expression of each constant c'p . In the two dimensional case, 
Le Gall [16] also showed that L{t) admits the asymptotic expansion in powers of 
1/logl 

When d is odd, Hamana [6] showed that the asymptotic expansion for L{t) 
with the help of fl4.2p . The purpose in this section is to improve the asymptotic 
behavior of L{t) if d is even and not less than 6. Throughout this section, we use 
Cj's for positive constants independent of the variable. 

Theorem 5.1. If d is even and not less than 6, there is a family of constants 
{a^n^Y^J^ such that 



' d-A 2'^/2-2(rf-4)r(c//2- l)t^/2-2 
1 '^a^^ r{{d-?>)/2)r^^-' logt I \ 

'2-1 tn/2 + _ 2)r{d/2 - 1)3 + ytd^z J 



td/2 

n=0 



where ^ d 2 
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Recalling the notation G'^'^^ = Gd/2-1, we set 



Li{t) 



L2{t) 



725 



,3 



e 2t 




-xy 



Vt Jo Jo y'G^'Ky) 



e 2* dxdy. 



Then we have from Theorem 14 . 1 1 1 hat . if d is even and not less than 6, 



L{t) = Sd-ir 



d-2 



{d-2)t ^ r2 



U{t) + {-lY'^-^L,{t) 



The calculation of Li{t) is easy since Re(2;j'^^) < for each j = 1,2 
X ^ and an integer n ^ we put 

RJx) = e~ 



(5.2) 
iVd. For 



E 

fc=0 



-xK 



and let M be a positive integer. Then it follows that 

M 







dx = 



n„2n /-oo 



n=0 
A/ 



n!(2t)« 7o 

(-l)V2"(2n)! 1 
^2-n!(zf))2"+it" 



1 



M 



2 2 > 



2t 



+ 



Hence we obtain 



2 ^ (-l)"r2"+3(2n)!C; 



2n+3 



n=0 



2"n! 



^n+1/2 



+ 



fM+3/2 



(5.3) 



as t — )■ 00, where the notation (jf^ is used to denote Cd/2-i,fc for an integer k ^ 1. 

For a proof of Theorem 15.11 we need to give an asymptotic behavior of L2{t) 
for large t. Let m = d/2 — 1 for simplicity and set -^^^(t) = L2(2r^t)/r^. We have 



LO(t) 



1 



00 /"oo e~^^ 



/Trtio Jo y'^Gm{y) 

We note that Gm{x) = Km{xY + n'^Im{x)'^. 
Lemma 5.2. For an integer n ^ 1 we have 



e 4t dxdy. 



(5.4) 



i?„_i(x^ 



2n. 



-c/x 



-l)"7r 



2r(n+ 1/2)' 



(5.5) 
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Proof. From a change of variables from x to y given hj y = x'^ we deduce that 
the left hand side of (15. 5p is equal to 



2 yn+l/2 ]^ k\ 



(cf. [al p.361]). By the formula r{z)r{l - z) = 7r/sin(7r2), ^ G C \ Z (cf. [H 
p.3]), we have 



^2 y r(n + 1/2)' 
which yields fl5.5|) . □ 

We give several constants which we need to describe the asymptotic behavior 
of 1/Gm{x) as X J, 0. For integers /i, k with l^/i<A;^m — Iwe put 



ki+k2A Vkh=k 



where 



We set 



^ r{m-h)r{m-k + h) 
k - /^kpi^^yl^^ r{h + l)r{k-h + l) ' 

k 



}2bk,k {k = l,2,...,m-l), am= 1 . 



We note that ai =bi. Moreover recall that = 1/4"^ ^r{m)'^. See fl2.25p for the 
definition of Hy. The second lemma gives the asymptotic behavior of 1/Gm{x)- 

Lemma 5.3. We have that, as x 10, 



G„,ix 



2 r m—1 



. m-l >. 

y + Yl + - + Oix^"") I . (5.6) 



Proof. By the series expression of the modified Bessel function 7^, we have that 
Im{xY is of order x'^^ as x J, 0. It is known that 

fc=o ■ ^ ^ 

2k+m 



+ (_!)-+! log- V i ( -Y 

^2^ A;!(/t + m)! ^2^ 



29 



where -0 is the logarithmic derivative of the gamma function (cf. [201 P-80]). This 
formula immediately yields 



r(m) /2\'"!^' {-l)T{m-k) ( 



X 



2 Va^y ^ r{k + i)r{m) \2 

i (m + 1) \2y X 

J, 0. Hence we have 



X 



k=0 

Therefore we conclude that, as x J, 



4 \x J r{k + l)r{m) \2 



2m\ 



GUx)=KUxy + Oix 

1 



1 r 1 1 

-^1-5^ hx'' - a^x^- log - + Oix'n . 

fc=l ^ 



We set G^{x) = Kmx'^"^Gm{x) for simplicity. It is sufficient to obtain the 
asymptotic behavior of 1/G^{x). We can easily derive 

^ m-l ,m-l 

In the case of m = 2, (15. 7p immediately implies (15. 6p . We concentrate on consid- 
ering the case of m ^ 3. In this case, the summation in the right hand side of 
dSZD is 



1 + $^(5^6^x2') +0^x2"^ log - + 0(x2"^) 



{5.t 



A simple calculation shows that the double sum in the right hand side of (15. 8p is 
equal to 

m— 1 771—1 

$^$^6.,.x^' + 0(x^™). 

h=l k=h 

Hence we deduce 

1 """^ 1 

—-- = l + Y. a,x'^ + a^x'- log - + 0(x^-), 

which implies (15. 6 p for m ^ 3. □ 
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We now proceed to a proof of Theorem 15.11 From (15.41) it follows that 

Q POO POO -| 

v^rjo Jo y^Gmiy) 



which is the sum of 

m-2 



O '"^ f ( 1 \n poo poo 1 

r\ POO POO 1 

We easily derive 

^ n=0 " 



4"n! t"+l/2 |/2n+4(^^^^) 



-y 



^ n=0 



where pj.''^ = ^^^fc. 

For X > and an integer /c = 0, 1, 2, . . . , m — 1 let 



1 



^ ' n=0 



k 

2n 



where we have put oq = 1 for convenience. We need the following lemma to derive 
the large time asymptotics of I^Kt). 



Lemma 5.4. We have that (t) is the sum of the following three integrals; 

2 roo r-oo ^-2\/tuy 

— — Kmy'^'^yaky'^^Rm+k-2{u'^)dudy, (5.9) 



k=0 

m— 1 



, r r (zM^M,„.„..-.,«,, (5,0) 

Jo Jo y ^ (m + A:- 1)! 

2 /-oo /-oo ^—2\/tuy 

-Qm-i{y)R2m-2{u^)dudy. (5.11) 




Proof. By the definition of Qo and -R^, we have 

2 /-oo /-oo g-2v^u{/ 



L^(t) =^ y y ^3 Kmy'^'^Rm~2{u^)dudy 



^ Jo Jo y^ (m- 1)! 



2 /-oo /.oo g-2v^uy 

+ ^ / / ^ — Qo{y)Rm-i{u^)dudy. 

\/'^ Jo Jo y 
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Moreover, it follows that, for y > 0, u > and an integer k ^1, 

/I \m+k-l 



{m + k- 1)! 

Taking the sum on k over [l,m — 1], we deduce that the third term of the right 
hand side of (15.1 21) is equal to 

2 roo POO ^-2Vtuy 
— / / § Q^rn-l{,y)R2m-2{u^)dudy 

Jo Jo y 

2 roo roo ^—2\/iuy m—l 




Jo Jo 



f^mV'^"' ^ aky'^^Rm+k^2{u'^)dudy 



k=l 



Hence we conclude that Ll{t) is the sum of fl5.9l) . flS.lOp and (15. lip . □ 

For (15. 9 p we first carry out the integral in y and use (15.51) . Then we see that 
(15. 9p is equal to 



2AC^'^4 afc(2m + 2fc-3)! Rrn+k-2^, 

^ 2-^ 22m+2k-2fm+k-l ^2(m+fc-l) 



fe=0 



-l)'"+'=-iafcr(2m + 2fc - 2) 1 
22™+2fc"2^(m + A; - 1/2) t^+^^i 



(5.13) 



Since 



r(2m + 2A; - 2) 2'^"'+'^^-^r{m + k - I) 



r{m + k- 1/2) 

which is the direct consequence of (14. Sp . the right hand side of (I5.13P and so (15.9 
are equal to 



~9~ 



fc=0 



2 ^ ^m+fc-l 
fc=0 

Carrying out the integral on u in (I5.10p . we see that (I5.10p is equal to 



2 !^^ (-i)rn+fc-i(2m + 2fc-2)! 1 Q^jy) 

0f2^ 22"^+2*^-l(m + fc-l)! t"^+fc-l/2 ^2m,+2fc+2 ^' 



which coincides with 

m—l 



1 (-l)-+^-^r(m^ + A:-l/2) /-^ gfc(y) 

^ 2-^ frn+k-l/2 ^2m,+2fc+2^' l^'-'-^J 



k=0 
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Here we have applied 

r{2z-i) 1 r{2z) 22^-2 / 1 



r(z) 2z - 1 r(z) 



TT 



riz--]. (5.15) 



We should note that the integral in the right hand side of fl5.14p converges for 
each integer k = 0,1,2, ... ,m — l. It is easy to see that fl4.13l) immediately yields 

that, as X — oo, 

Qk{x) = «:^x2-+2'= + o(x2™+2'=). (5.16) 
Moreover we deduce from (15.61) that, clS X 

ffi:^afc+ia:2'"+2'=+2 + o(a;2m+2fc+4) if ^ A; ^ m - 2, 

Qk{x) = < 1 

KmO-mX^^ log — h 0(x^"^) H k = m — 1. 

K X 

Hence Qk{y)/y'^"^^^''~^'^ is integrable on (0, oo) for each = 0, 1, 2, . . . , m — 1. 
For t > we let 

Ji y Jo 



oo 



Jo y 



Psit) = fy'-~Hog-dy r e-'^'-yR,^.,{u')du, 

where we put 

Qm{x) = Qm~i{x) - K^a^a^'^^log-. 

X 

Then flHTT]) is the sum of Pi(t), P2{t) and P^it). 

By virtue of (I5.16p . we obtain that \Qjn-i{y)\ ^ Ca?/''™"^ for y ^ 1. Combining 
this estimate with \R2m-2{u'^)\ = C^u'^"^~'^, we deduce 

\m\ ^ C, y^-^dy e-^^^-yu-^-^du = ^ f ^. 

This means that Pi{t) is of order i/i2m-i/2 jg negligible. 

We next show that P2{t) is of order l/t^™-!. It follows from (15. 6p that 
\Qmiy)\ = Cjy'^^ for y < 1. Noting R2m-2{x) ^ for x ^ 0, we obtain by 
the Fubini theorem that 

/•oo noo 

\P2{t)\ ^ Cs / -R2m-2iu^)du / c'^^'^y y^^^' dy . 
Jo Jo 

The formula (15. 5p implies that 

1^(^)1 ^ / -TZ—,, du 



t2™-i 7p ^^'"-a 2r(2m- l/2)t2m-i' 
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The calculation of Psit) is slightly complicated but not difficult. A change of 
variables from y to v given by 2\/ty = v yields 

For t > 1 /4 we set 

poo roo 

P^{t) = log(2v^) / v^'^-^dv / e-™i?2™-2(M')rfM, 
^0 Jo 

POO POO 

Pl{t) = log(2v^) / v^'^-^'dv / e-"''i?2m-2(n')dn, 
v^'^-^Xogvdv J e-^''R2m-2{u')du, 

Psit) = f log - dv f e-''^R2m-2{u^)du, 

Jo ^' Jo 

Jo ^' Jl 



Then we have 



We show that P^it) is the leading part of Psit) and the others are all negligible. 
Recall that R2m-2{u'^) = 0. It follows from fl5.5p and the Fubini theorem that 

Applying f l4.8p for z = 2m — 1, we have 

P.^{t) = -2^'""^v^r(2m - 1) logt + 0(1). 
The estimate of Pi{t) is easy. Indeed, we have 

dv 



POO POO POO 

|P|(t)| ^ Clologt / v^"'-'^dv e-^V^-'c/M ^ Ciilogt / 

J2Vt Jo J2- 



9 ' 

2Vt ^ 



which is of order \ogt/\/t. The way of estimates of the remaining integrals is 
similar to that of P^{t). We deduce 

\Pl{t)\^Ci2 v^'^-Hogvdv e-™u^™^2ciu ^ Ci3 / ^dv 
Jl Jo Jl ^ 

and that, by ^ g-™^4m-2 ^ ^ ^, ^ 

\PUt)\^Ciif v^'^-^log-dv ! e-^^M^'^-^du ^ Ci4 /" v^'^-'^log-dv. 

Jo V Jo Jo V 
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These immediately imply that -Pf (t) and P^it) are of order 1. Note that Rn{x) 
is asymptotically equal to (— l)""'"^x"/n! as x — )■ oo for an integer n ^ 0. This 
yields that |-R2m-2(w^)| is bounded by Ci^u^"^'^ for m ^ 1 and then we deduce 



\P^(t)\^Ci5 f v^^'^Hog-dv [ e-™M"cin ^ Ci6 / \og-dv = C 
Jo ^ ii Jo V 

Therefore, by virtue of Lemma 15.41 we accordingly obtain 
,2... _Km ^ (-l)"^+^-^r(m + k-l)ak 

k=0 

i!^^ (_i)^+fe-ir(m + fc-i/2) r_Qk{y)_, 

~^ ^ 2-^ ^m+fc-l/2 y2m+2k+2 y 

Kraamr{2m-l) \ogt ^ ^ / 1 



4 ^2m-l y ^2m~l 

which implies that we finished to give the asymptotic behavior of (t). 
Recall the definition of and am- We deduce from fl5.15p that 

Km.amr{2m-1) 1 (-1)" r(2m-l) 



4'"-ir(m)2 4'"-imr(m)2 



4"^i/7rmr(m)3 \^ 2^ 
Since L2(t) = r'^Ll{t){t/2r'^), we obtain, by ([53]) for M = - 4, that 

-^ (-l)"(2n)!r^"+3^£)^3 1 
7r ^ 2'^n! t"+V2 

^^=0 

2'^^^ (-l)'^/^+"(2n)!r^"+^gS+3 1 

n=0 

(-l)"r(rf/2 + n - 2)an2^/^+^-\^+^^-^d^'^ 

Hd/2-l ^d/2+n-2 
n=0 

r{{d-3)/2y'-' \ogt , ^/ 1 



where 



~Jo y^^^^ 



It follows from (El]) and (O that Q_i = and 
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for n = 0,l,2,...,(i/2 — 3. Therefore we can conclude that 
-Li(t) + (-l)'^/2"%(t)^ ^° 



1 ^ On 

~'~ A-,' ^n/2 

n=0 

^ r((rf-3)/2)r^'^-^ logt ^ ^/ 1 



Our proof of Theorem 15.11 is completed. 

6. Zeros of Macdonald functions 

We can find enough properties concerning zeros of J^, and ly (cf. [HI [20]). 
However there is less information on zeros of Ky. See [13] and [22], for example. 

Our purpose in this section is to represent all zeros of as the root of a 
polynomial of order N{v). Since Ky = K^^ and N{v) ^ 1 if \v\ ^ 3/2, it is 
sufficient to consider the case of ^ 3/2. Moreover, ifz/ = n + 1/2 for an integer 
n ^ 1, the formula 

yields that all zeros of Ky are the solutions of the equation 
Here we have used the notation 

r{u + k + 1/2) 



k\r{u-k + 1/2)' 



From now on, we discuss the cases when > 3/2 and v — 1/2 ^TL. By virtue 
of (12.61) . we have that, for x > 

, . = 1 H h > hcos(7rz/) / — ,^ , (6.1) 

We can derive the power sum of Zy^i, Zy^2, ■ ■ ■ , Zu,n{u) with the help of (16. ip . The 
Newton formula (cf. [H p. 276]) gives the polynomial whose roots are Zy^i, Zy^2, ■ ■ ■ , Zu^n{v)- 

Lemma 6.1. Let 

N 
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for positive integers n, N and zi, Z2, ■ ■ ■ , zn G C. We define a sequence {^n}n= 
of complex numbers by C,o = I and 



1 

= / ^ in-kPk 



n 

k=l 

for n = 1,2, . . . , N . Then we have that, for z E C 

N N 

n(^-^.)=E^^-^"- (6-2) 

j=l n=0 

Proof. Let s„ be the elementary symmetric polynomial of degree n, that is, 



lSjl<j2<-<jnSN 



Zj^ Zj^ . . . Zj^^ . 



The Newton formula yields that Pi = —si and 

n-l 

p^ = Y^{-lY-^+hn.kPk + {-iT^'ns^ 

k=l 

ioT n = 2,3, . . . , N. Therefore we easily deduce (6.2) from the formula 

N N 
j=l n=0 

This completes the proof of this lemma. □ 

We first consider the asymptotic expansion of (16.11) for large x. Recall that, 
for z/ ^ and a given integer M ^ 



* n=0 



(2x)" \x^+^ J 



as X ^ oo (cf. [H p. 123], [201 p.202]) and we have 



'^'^^o(-i^), (6.3) 



KJx) ^ \x 

^ ' n=0 ^ 

where {a„}*{l|)^ is the sequence of real numbers defined by 

(i^ + l,n) -^-^ {y,n — k) 



2n-fc 
A;=0 



for n = 0, 1,2, ...,M+1. A simple calculation shows 

1 1 / 2 1\ 1 A 9 1 

ao = l, ai = z/+-, a2 = -U-^l, «3 = -2(^^-4 

1^2 ^\( 2 25\ 1/2 l\/2 13 
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The remaining constants Oe, ay, . . . have comphcated forms. It is easy to give the 
asymptotic expansion of the right hand side of (16.11) and then we have that, for 
M ^ N{u) 

K+i{x) ^ 2u 1^ 1 cos(7rz/) f°° dy 



2iy 1 ^ 1 cosfyrz/) r°° 
1 + > ^ + ^ — - 

dj dj 1 ^y-:J I ^ ^ Jo 



and 

KJx) X ^ ^ 

A'/ ^ ^ 



Here we should note that the integral of y'^~^/Gp{y) over (0, oo) converges for 
each integer m ^ 0, which can be shown by Lemm [2l3] and fl2.25p . Comparing 
the corresponding coefficients in ( 16. 3 p and f l6.4p . we obtain 

iV(.)^.-i + cosMf (6.5) 



^ ^i:,^. = -a„+i + (-l)"cos(7ri.) / yrr-^dy (6-6) 

~[ Jo '-'iy[y) 



forn = 1,2,..., M. We define a sequence {an}n=o of complex numbers by = 1 
and 

1 " f /""^ y'''^^ 1 



for n = 1, 2, ... , N{v). Therefore, by (16. 6 p and Lemma I^TTl we have the following 
theorem. 

Theorem 6.2. For \v\ > 3/2 the zeros of are the solutions of 

n=0 

We obtain another polynomial whose roots are z^^i, z^^2, ■ ■ ■ , Zu,Niu) by consid- 
ering the asymptotic behavior of (16. ip for small x, which is an improvement of 
Lemma Ol Let ^ 1 be an integer with N + 1/2 < u < N + 3/2. It follows 
from Lemma [4.41 that, as x 4 0, 

^ ' n=0 
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where {bn}n=o ^ sequence of real numbers defined by 

u + l/2\2T{2iy -n + 2) fu - l/2\2''~^r{2u - n + k) . 



n J r(2z/ + 2) p^V^-^ J r{2u) 

for n = 0, 1, 2, . . . , 2N + 1. A simple calculation shows 

1 



6o = 1, 6i = 63 = 65 = 67 = 0, 62 



4z/(z/ - 1) 

^6 



16z/(z/- l)2(z/-2)' ' 32i/(j/- l)3(i/-2)(z/-3)' 
The remaining coefficients have complicated forms. We easily get 

2N N{u) 



[x] 2v 



KJx) X ' ^ + £^"1^ .n+l 



n=0 j=l ^'^J 



(6.8) 



We have to remark that we can not derive the higher term in f l6.8p while M in 
(16. 4p is arbitrary. This deference is caused from the integrabiliy of l/y'^Gj^{y). 
Comparing the corresponding coefficients in (16.71) and (16. 8p . we obtain 

g_.,.,„„_,r„.(..)^ _|_ ,.10, 

for n = 2, 3, ... , 2N. We define a sequence {/3^}^io'' of complex numbers by 

/3o = 1, 



/3'( = 1 + cos{niy) 



dy 



y^G^iy)' 

dy 



+-/3:-i|i+cos(H r 

^ I Jo 



y'G^iy) 

for 77, = 2, 3, ... , N{u). By (16. 9p . (I6.10p and Lemma 16. ![ we have the following 
theorem. 

Theorem 6.3. For \v\ > 3/2 the zeros of are the solutions of 

N{v) 
n=0 
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Remark 6.4. It is known that, if z/ — 1/2 is not an odd integer, 



I 71 



(cf. [201 p. 512]). Here 6^, is the unique number determined by 



COS 6', 



sin TTZ/, sin = cos vrz/. 



Hence we deduce from (16. 5p that, if z/ > 3/2 and z/ — 1/2 ^ Z, 

r dy ^ 9, 



(6.11) 



Jo yGu{y) 7rcos(7rz/)' 



Similarly, when ^ z/ < 3/2 and z/ 7^ 1/2, we can easily derive (16. lip by virtue 
of (12.51) and (16. 3p . This implies that, for an integer n ^ 0, 



which has been obtained in [22]. 

Remark 6.5. When 3/2 < v < 7/2, Niy) = 2 and the zeros of K^, satisfy some 
quadratic equations. In particular, when z/ = 2, the zeros 2:2,1 and 2:2,2 satisfy 



By using Mathematica we obtain from each system of equations that the zeros 
are close to —1.28 ± 0.43 z, and check the comment "The two zeros of K2{z) are 
not very far from the points — 1.29 ± 0.44 z. "in [20l p. 512]. Moreover we find that 
the zeros of are close to —1.68 ± 1.31 i in the similar way. 
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